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Abstract:  The  critical  hydraulic  configuration  for  a  set  of  torque-tube  gates  controlling  the  flow  through  the  navigable  portion  of  a 
spillway  was  experimentally  identified.  In  this  paper,  an  analytical  model  for  the  upstream  pressure  field  on  a  typical  gate  within  the  set 
is  constructed.  The  gate  rotation  from  the  maximum  elevation  (gate  in  closed  position)  and  the  hydraulic  torque  transmitted  by  the 
pressure  field  to  the  gate  tube  are  formulated.  Mean  values  of  parameters  of  response  are  often  sufficient  for  the  preliminary  design  of  a 
gate.  The  dispersions  of  these  parameters  of  response,  which  are  necessary  for  the  final  design  of  a  gate,  may  be  computed  using  the 
corresponding  mean-square  values.  These  were  obtained  empirically  in  a  flume  from  experiments  on  a  1/15-scale  physical  model  of  a  set 
of  three  prototype  gates  for  the  Montgomery  Point  Lock  and  Dam  project.  Theoretical  predictions  of  parameter  mean  and  mean-square 
values  compare  well  with  the  average  corresponding  statistics  obtained  experimentally. 
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Introduction 

The  flow  through  the  navigable  portion  of  a  river  barrier  may  be 
conveniently  controlled  by  a  parallel  array  of  torque-tube  gates, 
as  Fig.  1  illustrates  (see,  for  example,  Fletcher  and  de  Bejar  1995 
for  a  detailed  description  of  the  hydraulic  structural  10-gate 
scheme  in  the  Montgomery  Point  Lock  and  Dam  project). 

The  torque-tube  type  of  flap  gate  is  operated  by  controlling  the 
magnitude  of  the  torsional  moment  applied  to  the  structural  tube 
rigidly  connected  to  the  gate  panel  at  its  lower  horizontal  edge 
(Sehgal  et  al.  1997).  Fig.  2(a)  shows  a  typical  torque-tube  gate  for 
flow  control  in  navigational  spillways.  In  the  case  of  the  Mont¬ 
gomery  Point  Lock  and  Dam  project,  the  generally  curved  gate 
skin  adopts  a  circular  profile,  and  the  gate  elevation  is  identified 
by  the  angle  that  the  tangent  to  the  skin  circle  makes  with  the 
horizontal  plane  at  the  connection  to  the  tube.  The  gates  operate  at 
an  elevation  ranging  from  zero  (an  open,  horizontal,  position)  to 
70°  in  the  closed-gate  condition. 

Each  possible  combination  of  individual  operational  gate  el¬ 
evations  within  a  given  gate  set  is  subjected  to  complex  hydraulic 
loading,  and  the  task  of  identifying  the  critical  gate  configuration 
governing  the  hydraulic  structural  design  necessarily  involves  the 
application  of  both  analytical  and  experimental  techniques. 

Investigations  during  the  past  two  decades  have  directed  atten¬ 
tion  to  experimental  analysis  (Naudascher  1991;  Naudascher  and 
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Rockwell  1994)  and  to  updating  design  guidelines  of  spillway 
gates  (Sehgal  1996;  Sagar  1995).  In  particular,  Sehgal  (1996)  and 
Viparro  and  Hansen  (1993)  discussed  general  features  of  the 
torque-tube  type  of  flap  gate,  but  due  to  the  complexity  of  hydrau¬ 
lic  conditions,  recommendations  regarding  design  loads  are  lim¬ 
ited.  This  paper  presents  experimentally  verified  mathematical 
models  that  address  this  need  for  the  case  of  navigable  spillway 
torque-tube  gates. 

An  analytical  description  of  the  pressure  distribution  on  the 
upstream  face  of  a  typical  gate  in  the  critical  configuration  is 
derived,  and  the  associated  torque  transmitted  to  the  gate  tube  is 
computed  [Fig.  2(b)].  The  theoretical  pressures  and  torque  are 
compared  to  the  laboratory  data  from  experiments  on  a  carefully 
instrumented  1/15-scale  physical  model  of  a  prototype  Montgom¬ 
ery  Point  gate  at  the  U.S.  Army  Engineer  Research  and  Develop¬ 
ment  Center,  and  the  second-moment  characterization  of  the 
tube’s  torque  response  to  the  imposed  loads  is  verified. 

During  the  experiments,  it  was  recognized  that  the  largest  hy¬ 
draulic  forces  acting  on  a  single  central  gate  occur  when  all  gates 
in  the  experimental  flume  simultaneously  face  the  flow  in  a  fully 
closed  position  (de  Bejar  1995).  Fig.  3  illustrates  this  critical  con¬ 
figuration.  This  figure  shows  the  basic  depths  recorded  in  the  flow 
overpassing  the  system  of  gates:  the  depth  of  the  nappe  at  the  gate 
crest  (hf=  59  mm),  the  smallest  depth  in  the  free-falling  cascade 
(33  mm),  and  the  depth  of  the  enclosed  tailwater  pool  (128  mm). 
For  the  purpose  of  evaluating  design  hydraulic  forces  on  the  gate 
with  a  ventilated  nappe,  it  is  conservative  to  ignore  the  hydro¬ 
static  pressure  field  on  the  back  of  the  gate.  This  tailwater  pres¬ 
sure  was  verified  as  minor  for  other  gate  configurations  in  the 
experiments,  as,  for  example,  when  the  central  gate  in  a  set  of 
three  is  fully  closed  and  the  adjacent  lateral  gates  are  in  a  full 
open  position,  or  when  only  the  central  gate  in  a  set  of  three  is 
fully  open  and  the  adjacent  lateral  gates  are  in  a  full  closed  posi¬ 
tion. 

At  the  outset,  some  limitations  of  the  mathematical  models  to 
be  developed  are  noted:  (1)  The  change  in  the  flow  pressure  field 
due  to  the  gate  vibrations  is  neglected  and  (2)  the  models  accu- 
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Fig.  1.  Montgomery  Point  Lock  and  Dam:  application  of  torque-tube 
gates  for  flow  control  in  a  navigable  spillway 


rately  represent  the  behavior  of  the  set  of  gates  in  a  uniformly 
closed  position  (70°  elevation)  only  if  the  back  of  the  gates  is 
ventilated.  For  unventilated  gates,  strong  turbulence  may  form 
against  their  backs,  thus  invalidating  the  underlying  assumptions. 


Hydraulic  Pressure  Field 

To  construct  a  mathematical  model,  the  following  simplifying  as¬ 
sumptions  are  postulated.  Water  is  assumed  an  incompressible, 
inviscid,  and  homogeneous  fluid,  and  its  flow  is  taken  as  irrota- 
tional.  The  critical  hydrostructural  configuration  is  considered, 
i.e.,  the  upstream  flow  at  its  maximum  operational  headwater  pool 
running  over  the  torque-tube  gate  in  a  closed  position. 

Consider  the  idealized  unit-width  control  volume 
D' DA''' ABB' A" A'" D'  in  Fig.  2(b)  and  its  natural  partition  into  two 
subvolumes  by  the  vertical  plane  A'"An.  For  simplicity,  the  flow 
was  assumed  to  enter  and  exit  the  volume  with  uniform  velocities 
U  and  m,  respectively.  However,  this  assumption  for  the  exit  ve¬ 
locity  u  is  probably  unconservative  since  a  slightly  inclined  ve¬ 
locity  vector  generally  produces  a  larger  resisting  moment  at  the 
gate  tube.  The  forces  acting  on  the  boundary  surface  are  sche¬ 
matically  represented  in  the  diagram.  The  weights  of  water  in  the 
various  parts  of  the  control  volume  are  designated  by  Wh  i=  1,  2, 
and  3.  The  vertical  and  horizontal  components  of  the  force  ex¬ 
erted  by  the  gate  on  the  flow  act  at  point  Q  and  are  denoted  by 
and  R2,  with  lever  arms  r1  and  r2,  respectively,  about  the  pivoting 
point  A. 

Assuming  that  the  flow  through  AA"  is  approximately  horizon¬ 
tal,  the  vertical  equilibrium  of  control  volume  ABB' A" A  may  be 
expressed  as  (Prandtl  and  Tietjens  1957;  Shames  1995) 

R \  -W2-W3  =  ()  (1) 

Notice  that  the  discontinuity  A" A"’  in  this  geometric  simplification 
also  introduces  some  error.  In  reality,  the  flow  free  surface  pro¬ 
vides  a  smooth  transition  from  line  D'A'"  to  line  A"B' ,  as  indi¬ 
cated  schematically  by  the  dotted  curve  in  Fig.  2(b).  Therefore,  in 
reality,  there  is  a  small  additional  water  weight  acting  on  A"B' 
that  is  neglected  in  the  geometric  idealization. 

In  the  absence  of  horizontal  body  forces,  the  equation  of  con¬ 
servation  of  linear  momentum  applied  to  the  whole  control  vol¬ 
ume  in  Fig.  2(b)  may  be  expressed  as 


Fig.  2.  (a)  Typical  torque-tube  gate;  (b)  idealized  control  volume 


Pl-P2  —  R2  —  —  U2pAv  +  u~pAB  (2) 

where  Px  and  P2=resultant  forces  of  pressures  acting  on  sides 
DD'  and  AA1V  of  the  control  volume,  respectively;  AD  and  A B 
= cross-sectional  areas  of  the  control  volume  at  Locations  D  and 
B,  respectively;  U  and  u= corresponding  flow  velocities,  assumed 


Fig.  3.  Physical  model  of  Montgomery  Point  torque-tube  gate.  De¬ 
sign  controlling  condition:  lateral  view  of  three  gates  in  the  experi¬ 
mental  flume  in  closed  position 
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Fig.  4.  Area  subtended  by  gate  arc  AB  (shaded) 


to  have  a  direction  normal  to  the  respective  cross-sectional  areas; 
and  p=mass  density  of  water. 

Notice  that  the  assumption  of  a  zero  hydrostatic  pressure  on 
the  segment  BB'  implies  that  this  segment  is  taken  a  short  dis¬ 
tance  into  the  free  jet.  The  approximation  of  a  zero  gauge- 
pressure  free  jet  is  adequate  for  small  values  of  hf  (Street  et  al. 
1996). 

The  theorem  of  moment  of  momentum  about  point  A,  applied 
to  the  whole  control  volume,  produces  the  following  equation: 


H 


—  -  a  sin  60  +  WA— (cos  0O  -  cos  0y)  +  W3  ■  x3  +  P2  ■  y2 


-  Ri  ■  r3  -  R2  ■  r2  =  -  I  —  -  a  sin  0O  I  t/2pAr 


H 

2 


hf~ 

a(stn  0y- sin  0O)  +  itpAB 


(3) 


where  H  and  approach  depth  and  the  brink  depth,  respec¬ 
tively;  0O  and  0y=elevation  angles  of  points  A  and  B,  which  are 
the  respective  ends  of  the  circular  gate  AB  in  closed  position;  and 
a = gate  radius. 

Since  point  Q  is  on  the  circle,  the  following  relationship  is 
satisfied 


( a  cos  0O  -  tq)2  +  ( a  sin  0O  +  r2)2  =  a2  (4) 

Next,  consider  the  generating  radius  vector  with  center  at  C,  mag¬ 
nitude  a,  and  argument  0  e  (0o,0d,  shown  in  Fig.  4.  The  shaded 
area  A^  is  given  by 


A  = 


ra  cos  Ca  sm  “f 

_  dydx 

J  a  cos  0y  J  ytp-x2 

"2  sin  Qf( cos  0q  —  cos  0y)  +  “(sin  20y-  sin  20o) 


9/  ~  8q 
2 


(5a) 


and  its  centroidal  point  has  abscissa 


x  = 


a 3  sin2  0o(sin  0y-  sin  0O) 

(5b) 


The  weight  of  water  corresponding  to  the  area  As  is  W^y-Aj, 
where  y  =  p-g=the  unit  weight  of  water  and  g=  acceleration  of 
gravity.  The  line  of  action  of  W3  is  defined  by 


x2  =  a  cos  0O  -  x  (5c) 

Now,  the  reactive  force  R{  may  be  computed  from  Eq.  (1)  as 

R\  =  W2  +  W3  =  y  ■  a  ■  hf-  (cos  0O  -  cos  0y)  +  W3  (6) 

in  which  the  brink  depth  (hf)  may  be  approximated  as  that  at  a 
free  overfall,  and  this  may  be  expressed  as  the  coefficient  of  con¬ 
traction  (Cc)  times  the  hydrostatic  critical  depth  (hc) 

hf=Cc-hc=Cc-(^j  (7) 


where  q = rate  of  discharge  per  unit  width  of  gate  (Rouse  1936). 

On  the  other  hand,  the  resultant  forces  P,(i=  1 ,2)  can  be  esti¬ 
mated  with  the  expressions 


P  i  = 


y  ■  H2 

2 


(8a) 


and 


P2  = 


y  H2 
2 


1 


1  -  —sin  0n 
H 


m 


corresponding  to  the  hydrostatic  pressure  distributions  on  DD' 
and  AA‘V,  respectively.  The  approach  depth  (H)  may  be  obtained 
from  the  discharge  rating  curve  (e.g.,  Bradley  1954). 

The  flow  velocities  U  and  u  may  be  estimated  directly  on  the 
basis  of  continuity,  given  the  incoming  flow  q,  i.e. 


u=U  ■ 


q_ 

hf 


(9) 


Now,  the  reactive  force  R2  may  also  be  computed,  from  Eq.  (2), 
as 


R2  =  (Pl-P2)  +  1(U2-H-u2-hf)  (10) 

S 


Fig.  5  shows  the  free-body  diagram  of  the  idealized  gate  sub¬ 
jected  to  the  hydraulic  forces  Ri  and  R2  (per  unit  width  of  gate) 
and  the  corresponding  internal  reactions  at  the  base  of  the  can¬ 
tilever,  where  N0,  V0,  and  M0  represent  the  axial  load,  shear 
force,  and  bending  moment  per  unit  width  of  gate.  The  equilib¬ 
rium  formulation  leads  to  the  following  force  transformation 
(Przemieniecki  1985): 


cos  0O 

sin  0O 

M 

-  sin  0O 

cos  0O 

1  '<■  j 

(11a) 


which  allows  the  direct  evaluation  of  the  shear  force 


V0  =  Ri  sin  0o  +  R2  cos  So  (11  b) 


Eq.  (life)  provides  one  of  the  boundary  conditions  for  the 
boundary-value  problem  to  be  formulated  next.  The  other  bound¬ 
ary  condition  is  given  by  the  fact  that  the  curved  cantilever  has  a 
free  end  at  B.  Mathematically  expressed,  one  has 
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Vocos(Bo) 

Fig.  5.  Free-body  diagram  of  gate  section 


v(e0)  =  Vq 

V(Qf)  =  0 


(12) 


Now,  if  the  analysis  leading  to  Eqs.  (6)  and  (10)  for  the  hydraulic 
forces  is  generalized  for  an  end  elevation  0,  instead  of  0,  [Fig. 
2(b)],  the  following  expressions  are  obtained: 


7^(0)  =  7]  ah( cos  0O  -  cos  0)  +  a ■ 


cos  0o(sin  0  -  sin  0O) 


-  -(sin  20  -  sin  20o)  -  -(0  -  0O) 


(13a) 


and 


*2(e)  = 


2  U2I  1  1 


1  -  —sin  0n  -  2  ,  . 

H  7  g\h  Hi  J 


(13  b) 


where  h= water  depth  of  the  nappe  at  the  crest  of  the  generalized 
cantilever  AF.  After  some  algebraic  manipulations,  the  first  two 
derivatives  of  each  of  these  expressions  may  be  obtained  as 


dR  i  , 

- =  ya  cos  0 

dQ 


h  \  T|2(cos  0o  -  cos  0) 

-  tan  0  + - « - 

aj  r|2  -  1 


(14a) 


d2R  i  2 

- —  =  'vr 

dQ 


sin  20 

2  =7A“  -  COSQ  +  —y - 


2(ti2-1) 


_n 

ti2-1 

a  \  /  1  +  cos  20 


+  — — -(sin  20  -  sin  0  •  cos  0O) 


(cos  0O  -  cos  0) 


+7 


(V-i)3. 


(14b) 


dR2  T| 

— —  =  yah  cos  0— 


dQ 


T0  “I 


(14c) 


and 


Fig.  6.  Free-body  diagram  of  gate  differential  element 


dlR2_  J 

cos2  0  / 

dQ2 

[t]2-1 

V-i  \ 

sin  0 


3(ti  cos  0)2} 

(ti2-!)3  j 


(14  d) 


where 


ri2  = 


u- 


gh  I  \  h  I 


HV 


(14c) 


Next,  consider  the  free-body  diagram  of  the  differential  element 
of  gate  in  Fig.  6.  The  equilibrium  of  forces  in  the  radial  and 
tangential  directions  leads  respectively  to 


,  ,  dV 

R,  sin  0  +R7  cos  0  +  N=  — 
1  2  dQ 


(15  a) 


and 


,  ,  dN 

R2  sin  0  -  R,  cos  0  +  V  +  —  =  0 
2  1  dQ 


(15  b) 


where  the  prime  symbol  represents  derivative  with  respect  to  0. 

Eliminating  dN/dQ  from  Eqs.  (15a)  and  (15b),  one  obtains  the 
following  governing  second-order  ordinary  differential  equation: 


(15c) 


which  together  with  the  boundary  conditions  (12)  constitute  a 
boundary-value  problem  for  7(0). 

This  boundary-value  problem  is  iteratively  transformed  into  an 
initial-condition  problem  in  the  0  domain  and  subsequently 
solved  numerically  using  a  Fourth-Order  Runge-Kutta  method 
(Kreyszig  1988).  The  total  pressure  field  is  then  given  by  (see 
Appendix  for  the  derivation  of  a  closed-form  analytical  expres¬ 
sion  of  the  hydraulic  pressure  field  on  the  upstream  gate  skin) 


PrW  = 


1  dV 
a  dQ 


(16) 


The  calculation  of  the  torque  per  unit  width  of  gate  transmitted  by 
the  hydraulic  pressure  field  to  the  gate  tube  ( MR )  requires  the 
previous  evaluation  of  the  lever  arms  r1  and  r2.  These  lever  arms 
are  obtained  from  the  simultaneous  solution  of  Eqs.  (3)  and  (4). 
Segment  y2,  required  to  completely  define  this  system  of  equa¬ 
tions,  may  be  computed  from 
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forces 


Fig.  8.  Structural  coordinates  and  single-degree-of-freedom  (1-DOF) 

6i 


where 


1  y/L3 
y~P,  6 


V2  =  a  sin  0O  -  y  (17a) 

a  Y  ( a\(  a  V 
--stne0]  -31  —  Ml  -  —sin  0O  I  sin  0O 

(m 


the  centroidal  point.  The  torque-tube  reduces  to  point  A.  The  con¬ 
nection  of  the  operating  arms  to  the  base  of  the  gate  at  point  A  is 
rigid. 

The  displacements  are  discretized  in  this  model  by  the  three 
coordinates  shown  in  Fig.  8,  i.e.,  rotations  0,  and  02  at  Nodes  A 
and  D,  respectively,  and  the  linear  displacement  A  of  the  end  D  of 
the  operating  arms.  Only  the  strain  energy  associated  with  bend¬ 
ing  deformations  of  the  operating  arms  is  considered.  The  stiff¬ 
ness  matrix  of  the  system  may  be  formulated  as 


Subsequently,  the  hydraulic  torque  per  unit  width  of  gate  may  be 
obtained  as 


Mr  —  R\rl  + 


(18) 


Structural  Response 
Mean  Gate  Rotation 

The  gate  system  is  conceived  as  a  rigid-body  torque-tube  circular 
gate  standing  on  simple  supports  and  being  acted  upon  by  the 
hydraulic  pressure  field  and  by  the  viscoelastic  system  of  resisting 
forces  provided  by  the  operating  arms  and  the  horizontal  actuator. 

Consider  the  rigid-body  gate  in  closed  position  (elevation 
=2a)  shown  in  Fig.  7,  standing  on  simple  supports  and  being 
excited  by  an  external  dynamic  pressure  field  p  and  its  own 
weight  W  (mass  M=WI g).  The  overturning  moment  of  the  exter¬ 
nal  forces  about  the  center  of  the  tube  is  resisted  by  the  eccentric 
force  exerted  by  the  horizontal  actuator  with  stiffness  Ka  and 
equivalent  viscous  damping  Ca,  acting  at  the  end  of  the  operating 
arms  with  flexural  rigidity  El  and  length  L. 

The  analytical  model  for  this  structure  is  shown  in  Fig.  8.  The 
gate  is  idealized  as  a  rigid  circular  arc  (radius  a)  spanning  the 
elevation  0  from  0O  (at  A)  to  0^  (at  B),  and  its  mass  is  lumped  at 


2 


[K]  = 


2EI 

L 


1 


3 

L 


3 

1 

L 

3 

2 

L 

3  /  6  Ka  ■  L  ■  cos2  a  \ 
L  \L2+  2EI  ) 


(19) 


Elimination  of  coordinates  02  and  A  from  the  corresponding  stiff¬ 
ness  formulation  by  static  condensation  (Przemieniecki  1985) 
produces  the  generalized  stiffness  coefficient  for  the  degree  of 
freedom  0!  only  as 


K*  = - - -  (20) 

LI 

3 El  +  Ka-  L2 ■  cos2  a 

A  similar  procedure  applied  to  the  viscous  damping  matrix  ren¬ 
ders  the  generalized  viscous  damping  coefficient  associated  with 
the  degree  of  freedom  0[  only  as 

C*  =  Ca  ■  L2  ■  cos2  a  +  C*D  (21) 

where  C^= generalized  viscous  damping  coefficient  associated 
with  that  of  an  equivalent  viscous  damper  representing  the  energy 
dissipation  mechanism  from  the  structure  itself  plus  the  effect  of 
the  drag  force  on  the  gate  opposing  its  velocity  through  water. 
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Fig.  9.  Dynamic  equilibrium  of  1-DOF  rigid-body  gate  system;  vir¬ 
tual  displacement 


Fig.  9  shows  the  generalized  single  degree-of-freedom  system 
(cp  =  Oi)  vibrating  under  the  action  of  the  vertical  and  horizontal 
components  of  the  mean  hydraulic  force  on  the  gate  ( Rl  and  R2, 
respectively),  the  gate’s  own  weight  (VT)  and  the  resisting  elastic, 
viscous,  and  inertial  moments  (Ms,  Mc,  and  M,-,  respectively)  at 
the  gate  base.  Notice  in  the  figure  the  viscous  damper  with  coef¬ 
ficient  CD  that  exerts  the  drag  force  on  the  gate.  The  degree  of 
freedom  ip  is  defined  as  the  rotation  of  the  gate  with  respect  to  its 
static  weight-deflected  configuration  (Clough  and  Penzien  1993). 

The  dynamic  equilibrium  of  the  rigid-body  gate  is  formulated 
by  applying  the  virtual  work  principle  on  the  structural  system 
subjected  to  the  virtual  displacement  8cp  (i.e.,  to  a  small  variation 
from  the  instantaneous  position  of  equilibrium  defined  by  cp) 
(Mura  and  Koya  1992). 

Ignoring  negligible  terms,  the  energy  equation  2(5Work)=0 
reduces  to 


d2  cp  dcp  2 


Ryrx+R2r2  M0 

=  —  (22  a) 


h  h 

where  70=mass  moment  of  inertia  of  the  gate  about  the  tube  axis 


co 


2 


3  El 

l^L 


(22  b) 


the  square  of  the  natural  circular  frequency  of  vibration;  and 
=  C*/I0= overall  equivalent  viscous  damping  ratio  of  the 
upstream-wet  gate  system. 

After  taking  expectations,  the  steady-state  part  of  the  solution 
of  Eq.  (22 a)  (i.e.,  the  mean  steady-state  gate  rotation)  tends  as¬ 
ymptotically  to  the  limiting  value  (Clough  and  Penzien  1993) 


9l  = 


M  o 

cn2-/o 


where  M0=E(MR)= mean  value  of  MR. 


(23) 


Second  Moment  of  Hydraulic  Torque  on  the  Gate  Tube 

Let  the  hydraulic  pressure  field  be  monitored  by  nine  pressure 
cells  located  at  discrete  points  (strategically  distributed)  on  the 
upstream  face  of  the  gate.  The  governing  equation  of  motion  for 
the  gate  subjected  to  the  upstream  pressure  field  may  be  written  as 


d2v  i 
dt 2 


+  2U2tt/i) 


dvy 

dt 


(2tt/i)2Ui 


9 

■  2  (Ai  ■  fiiPiit) 

i=l 


(24a) 


where  Pj(i=  1 ,2, . . .  ,9)=pressure  acting  at  the  ith  pressure  cell 
location,  assumed  as  uniformly  distributed  over  its  tributary  area 
(A,);  fy  = natural  frequency  of  vibration  of  the  dry  gate  system; 
Vy=dytp  and  y^dj/dy,  with  c/,=distance  from  the  line  of  action 
of  the  force  associated  with  p;  to  the  tube  axis. 

Consider  the  particular  case  when  p,  =  8(f)  (i.e.,  a  Dirac  delta 
function  applied  at  the  starting  time)  acts  alone.  Eq.  (24a)  then 
becomes  (Bendat  and  Piersol  1993) 


d2vu  ,  .dvu  ,  ( dy\  ,  „  . . 

—T  +  2£w(2'n/i)—  +  (2'rrfl)2vu=  (y  j  •  (A;7i)8(f) 

(24  b) 

By  Fourier  transforming  Eq.  (24 b),  one  obtains  the  complex  fre¬ 
quency  response  function  between  the  pressure  p ,  and  the  dis¬ 
placement  vy  as 


Hu(f) 


2 


+  «’( 2Q 


(25) 


where  K’  =  3El/L. 

Therefore,  the  Fourier  transform  of  Vy  is  given  by  (Crandall 
and  Mark  1963) 


(26) 

where  P^f)  =  Fourier  transform  of  Pj(t). 

The  Fourier  transform  of  the  total  response  Vy(t)  is  conve¬ 
niently  obtained  by  superposition  in  the  frequency  domain  as 


9 

Vy(f)  =  2  Hu(f)  •  P,(f)  (27) 

1=1 


By  taking  the  inverse  Fourier  transform  of  Eq.  (27),  the  total 
response  Vi(t)  in  the  time  domain  is  obtained  as 


9  j2 

»,«>  -sS 


j2  r+co 

- 

i=  1  ^  J  -  oo  j 


P,(j)el^f,df 


-{7J  +Mly 


(28) 


The  total  hydraulic  torque  on  the  gate  tube  is  given  by 


MR(t)  =  K\(t) 


(29  a) 


or 


MR(t)  -  dy  ■ 


9  r+cc 

2  (^;7;) 

1=1  2  _oo 


1  - 


P  i(f)ei2™f‘df 


(29  b) 


If  the  pressure  pt  is  assumed  a  zero-mean  stationary  (and  ergodic) 
process  (by  extracting  the  time-average  from  the  record),  then,  on 
the  basis  of  linearity,  the  response  v , ,  and  its  contribution  to  the 
hydraulic  torque  (MR  l)= zero-mean  stationary  (and  ergodic)  pro¬ 
cesses.  Consequently,  the  following  relations  among  their  corre¬ 
sponding  spectral  density  functions  hold 
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=  1  j  )  -\Hu(f)\2-Spi(f) 


7  \  2 
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$»(/) 


1  - 


(2U 


/'  2 

fi 


(30) 


By  adding  the  contribution  of  all  p,’s  to  the  total  power,  one 
obtains  for  the  mean  square  of  the  hydraulic  torque  on  the  gate 
tube 


9  r 

E(M2R)  =  dlJJ(Aiyl)2- 

1=1  J -c 


i-i' 


(2U1  j 


(31a) 


9  /»+c 

E(M2R)  =  d2]-JJ(Aly,)2- 

1=1  Jo 


GM)df 


i-i' 


(2U(  j 


(31  b) 


in  terms  of  the  half-range  spectral  density  function  Gpi(f)  (Bendat 
and  Piersol  1993),  and  the  corresponding  variance  is  given  by 


aiK  =  E(M2r)  -  [E(MR)f  =  E(Mr)  -Ml  (31c) 


Experimental  Verification 

It  should  be  emphasized  that  multiple  flow-induced  loading  con¬ 
ditions  on  the  physical  model  were  examined  in  the  overall  study, 
corresponding  to  diverse  gate  operating  modes  (Fletcher  and  de 
Bejar  1995).  However,  the  theoretical  models  developed  in  this 
paper  consider  only  the  controlling  gate-system  configuration, 
whose  geometric  characteristics  were  idealized  for  mathematical 
tractability,  in  the  light  of  the  conservation  principles  of  fluid 
mechanics  and  the  theory  of  stationary  random  processes  for 
structural  dynamics.  Other  gate  and  flow  configurations  require 
different  formulations,  since  they  constitute  substantially  different 
physical  problems.  Among  the  complicating  factors  observed  in 
other  gate  configurations  are:  the  very  influential  effect  of  a  tur¬ 
bulent  pressure  field  on  the  downstream  face  of  the  gates,  exerted 
by  a  returning  tailwater  jet,  the  vibrations  induced  by  the  en¬ 
trapped  air  in  unventilated  nappes,  and  the  edge  effects  at  neigh¬ 
boring  partially  open  gates. 

Mean  Pressure  Field 

The  theoretical  formulations  are  verified  with  an  instrumented 
1/15-scale  brass  model  of  a  Montgomery  Point  torque-tube  gate 
(Fig.  3).  Fig.  10  shows  a  schematic  perspective  of  the  physical 
model,  indicating  the  location  of  nine  remotely  controlled  (ther¬ 
mally  insulated)  pressure  cells  used  to  monitor  the  hydraulic  pres¬ 
sure  field. 


Fig.  10.  Schematic  perspective  of  gate  physical  model 


The  geometry  of  the  gate  and  its  operational  hydraulic  envi¬ 
ronment  is  defined  by  the  radius  a =498  mm,  the  end  elevations 
0o=2O°  and  0^=72.75°,  the  depth  of  the  incoming  flow  d 
=  628  mm,  the  approach  depth  H=  557  mm,  and  the  design  rate 
of  discharge  per  unit  width  of  gate  q= 8  X  104  mm3/(s  mm).  Tak¬ 
ing  the  coefficient  of  contraction  as  Cc=0.68  [an  average  of  sev¬ 
eral  values  reported  in  the  literature  (Rouse  1936)],  Eq.  (7)  gives 
hf=  59  mm. 

The  mean  values  of  the  entering  and  exiting  flow  velocities 
may  be  calculated  using  Eq.  (9)  as  {7=0.144  m/s  and  u 
=  1.356  m/s.  Now,  Eqs.  (6)  and  (10)  give  the  vertical  and  hori¬ 
zontal  components  of  the  mean  force  exerted  by  the  flow  on  the 
gate  as  /?! =51 .28  N  and  R2=63.54  N,  respectively,  thus  allow¬ 
ing  the  calculation  of  the  boundary  shear  force  V0.  Eq.  (11/?)  gives 
F0=77.00  N. 

Given  the  numerical  solution  of  Eq.  (15c),  the  evaluation  of 
Eq.  (16)  yields  the  mean  pressure  field,  which  is  cast  in  nondi- 
mensional  form  by  normalizing  with  respect  to  the  approach  hy¬ 
drostatic  pressure  pn=y-H.  The  resulting  nondimensional 
pressure  distribution  (p,lp„ )  is  shown  graphically  in  Fig.  11  as  a 
function  of  the  normalized  circumferential  coordinate  s/La  on  the 
gate  arc,  with  origin  at  the  gate  crest  and  increasing  downwards 
(normalized  with  respect  to  the  total  arc  length,  Lfl=a-(0f-0o) 
=459  mm).  This  figure  also  shows  the  minimum  and  maximum 
values  of  the  experimental  pressure  acting  on  the  upstream  side  of 
the  1/15-scale  physical  model  of  the  gate  during  the  controlling 
100-s  test.  The  pressure  field  is  described  along  three  different 
gate  arcs:  the  left  edge,  the  middle  arc,  and  the  right  edge  (see 
Fig.  10).  Each  of  these  curves  represents  a  parabolic  fit  on  three 
data  points.  These  minimum  and  maximum  curves  delimit  narrow 
bands,  suggesting  very  little  fluctuation  of  the  pressure  at  a  given 
location  during  the  test.  The  theoretical  estimation  of  the  pressure 
field  is  an  excellent  upper  bound  over  most  of  the  domain  for  the 
gate  central  arc,  except  on  a  small  region  in  the  middle  third, 
where  the  estimation  is  slightly  unconservative.  Imperfect  sealing 
at  the  edges  of  the  instrumented  gate  (the  central  gate  in  the  set  of 
3)  allows  the  formation  of  lateral  jetting,  which  induces  an  ob¬ 
served  increase  in  pressure  on  the  upper  corners  of  the  gate,  as 
compared  with  the  theoretical  prediction  (left  and  right  edges  in 
Fig.  11). 
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Fig.  11.  Comparison  of  theoretical  prediction  of  hydraulic  pressure 
field  with  corresponding  experimental  pressure  range 


Mean  Torque  on  the  Gate  Tube 

The  numerical  solution  of  Eqs.  (3)  and  (4)  and  the  subsequent 
evaluation  of  Eq.  (18)  yield  the  mean  hydraulic  torque  on  the  gate 
tube  as  M0=  117.2  N  m/m.  This  theoretical  estimation  differs 
from  the  average  experimental  hydraulic  torque  during  ventilated- 
nappe  tests  (105.7  N  m/m)  by  11%,  on  the  conservative  side.  The 
discrepancy  between  the  theoretical  prediction  and  the  experi¬ 
mental  value  is  due  in  part  to  the  fact  that  the  quantity  experi¬ 
mentally  measured  is  the  bending  moment  at  the  end  of  the 
operating  arms  (sections  J  in  Fig.  10),  which  differs  somewhat 
from  the  hydraulic  torque  on  the  gate  tube.  The  agreement  be¬ 
tween  theory  and  experiments  is  generally  good. 

Mean  Gate  Rotation 

Fig.  12  shows  a  relevant  portion  of  the  free-vibration  decay  of 
the  physical  model  with  water  on  its  upstream  side  only,  and 
the  application  of  the  logarithmic  decrement  method  to  identify 
the  equivalent  viscous  damping  ratio:  £,w=9.5%,  on  the  average 
(Clough  and  Penzien  1993). 

If  Young’s  modulus  of  elasticity  for  brass  is  taken  as 
E=  106.9  X  103  MPa,  the  moment  of  inertia  of  both  operating 


TIME  [SEC] 


Fig.  12.  Identification  of  equivalent  viscous  damping  ratio  of  gate 
with  submerged  upstream  side 
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Fig.  13.  Spectral  density  of  bending  moment  at  cross  sections  J  of 
operating  arms,  i.e.,  at  the  joint  with  gate  tube 


arms  about  their  neutral  axis  in  bending  is  calculated  as  7=9.25 
X  104  mm4,  and  it  is  recognized  that  L= 236  mm,  then  Eqs.  (22 b) 
and  (23)  give  cpL=0.0315.  Experimental  measurements  on  the 
physical  model  with  a  ventilated  nappe  subjected  to  carefully  con¬ 
trolled  hydraulic  pools  gave  <pL= 0.0350°,  on  the  average.  There¬ 
fore,  the  theoretical  prediction  is  generally  in  good  agreement 
with  the  experimental  observations. 

Mean  Square  Torque  on  the  Gate  Tube 

Eq.  (29 b)  is  evaluated  numerically  using  the  pressure-cell  records. 
For  example,  using  the  controlling  test  with  a  duration  of  100  s 
and  data  obtained  at  a  rate  of  500  samples/s,  the  value  of  the 
power  adds  up  to  7s(M^)=4,054  (N  m)2,  which  is  compared  to 
the  experimental  value 

E(M2r)  =  (j2Mr  +  \i?MR  =  165.00  +  69. 172  =  4949  (N  m)2 

where  cr^= variance  of  the  random  process  MR ,  given  by  the 

area  under  the  corresponding  spectral  density  function  (Fig.  13), 

and  p,M  is  the  mean  of  the  process.  Again,  the  mean  square  of  the 
i  ^  /  *)  \  •  . 
directly  measured  process  E(MR)  differs  from  its  value  derived 

from  the  hydraulic  pressure  field  partly  because  the  directly  mea¬ 
sured  process  is  the  bending  moment  at  the  end  of  the  operating 
arms.  However,  results  compare  generally  well  (~18%  relative 
difference,  on  the  unconservative  side). 


Conclusions 

The  upstream  pressure  field  on  a  typical  navigable  spillway 
torque-tube  gate  is  constructed  analytically.  The  critical  hydro- 
structural  configuration  is  considered,  i.e.,  the  upstream  flow  at  its 
maximum  operational  headwater  pool  running  over  a  set  of 
torque-tube  gates  in  a  fully  closed  position.  Experimental  obser¬ 
vations  indicate  that  this  critical  condition  is  conservatively  sim¬ 
plified  in  the  theoretical  model  by  neglecting  the  effect  of  the 
tailwater  pressure.  The  mean  steady-state  rotation  response  of  the 
gate  is  formulated,  and  the  second-moment  characterization  of  the 
torque  transmitted  to  the  gate  tube  is  presented. 

The  mean  pressure  field  gives  an  excellent  upper  bound  over 
most  of  the  domain  for  the  central  region  of  the  gate  skin.  On 
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the  gate  corners,  however,  the  theory  underestimates  the  obser¬ 
vations.  At  least  in  part,  this  underestimation  is  due  to  edge  ef¬ 
fects  in  the  subscale  physical  model.  Derived  mean  values  of 
parameters  of  response  are  sufficiently  close  to  the  experimental 
observations  (on  the  average,  11%  overprediction  for  the  hy¬ 
draulic  torque  induced  in  the  gate  tube)  and  may  be  used  in 
the  preliminary  design  of  a  gate  for  the  application  being  consid¬ 
ered. 

The  dispersion  of  parameters  of  response  is  needed  for  the 
final  design  of  a  gate.  These  statistics  are  directly  related  to  the 
corresponding  mean  square  quantities.  In  this  study,  the  mean 
square  value  of  the  parameter  under  consideration  is  obtained 
empirically  (e.g.,  using  physical  modeling).  The  mean  square 
value  of  the  hydraulic  torque  induced  in  the  gate  tube  is  under¬ 
predicted  (on  the  average,  with  18%  error)  when  derived  from 
pressure  records,  as  compared  to  the  corresponding  mean  square 
value  measured  directly.  The  simplifying  assumptions  in  the 
structural  model  and  measurement  errors  are  the  major  contribu¬ 
tors  to  this  discrepancy.  In  practice,  however,  having  only  the 
pressure  records,  the  analyst  may  derive  straightforwardly  the  dis¬ 
persion  of  any  other  parameter  of  interest  without  having  to  mea¬ 
sure  it  directly. 

The  agreement  between  theory  and  experiments  is  remarkable, 
given  the  errors  associated  with  the  simplifying  assumptions 
under  the  theory,  and  the  design  of  the  physical  model  and  its 
actual  construction,  assemblage,  and  installation. 

The  physical  subscale  replica  of  the  scenario  for  the  torque- 
tube  gate  in  service  constitutes  an  imperfect  Froude’s  model  that 
allows  the  estimation  of  the  basic  design  hydrodynamic  forces  on 
the  prototype  gate  on  the  basis  of  similarity  principles  (Fletcher 
and  de  Bejar  1995).  Although  some  scale  effects  are  ignored  (like 
the  added-mass  effects  in  the  hydrodynamic  response  of  the  pro¬ 
totype,  for  example),  the  results  provide  invaluable  quantitative 
guidance  to  the  designer  that  would  be  difficult  to  estimate  with¬ 
out  the  experimental  observations  and  the  subsequent  analysis 
conducted  in  this  investigation. 
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Appendix.  Hydraulic  Pressure  Field 
on  the  Upstream  Gate  Skin  in  Closed  Form 

Governing  differential  Eq.  (15c)  for  the  shear  force  function  U(0) 
is  given  by 


V'+V=(2-R[+  R'2)  •  c0  -  (2  ■  R'2  -  R")  ■  50  (32) 

in  which  c0=cos(0)  and  50  =  sin(0).  Inserting  Eqs.  ( \Aa)-{\Ad ) 
into  Eq.  (32)  renders 


V"  +  V=ya2-k- 


3  h  1 

—  —  —  ^ ■ ^20  + 
2  a  r|2 


k2c0 


X  1  +  - 


2t0 


2  •  (tT  -  2)  +  — j 

2t0  . 


■  c20 


C0O  ,  .  3 k2  , 

+  —  •  (1  +  3  ■  c20)  + - •  50  ■  (c0o  -  c0)  •  c20 

2  h 

T)  ■  - 


(33) 

in  which  k=  1/(1  -  1/t|2).  The  general  solution  of  Eq.  (33)  is  ob¬ 
tained  as 


V(0)  =  A  ■  c0  +  B  ■  50  +  Vp(6)  (34a) 

where  the  particular  solution  is  given  by 


c0o  + 


^■0 
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■5  20 
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4r)2 
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(34  b) 


Upon  imposing  the  boundary  conditions  stipulated  by  Eqs.  (12), 
the  set  of  constants  {A,B}  in  the  complementary  solution  to 
Eq.  (33)  is  obtained  as 


M  i  |  [yo-v/0o)]- 50^X0/) -50o  | 

{B  j  “  sin(0/-  0O)  '  [-  [V0  -  VP(0O)]  ■  cdf-  Vp(6f)  •  c0o  J 

(35) 

Thereafter,  Eq.  (16)  provides  an  analytical  expression  in  closed 
form  for  the  total  hydraulic  pressure  field  on  the  upstream  gate 
skin 


1  dV(Q)  1  ya-k 
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(36) 
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When  the  approach  velocity  U  is  very  large  (U— >°°),  Eq.  (33)  for 
the  shear  force  function  V(0)  reduces  to 

V  +  y  =  _  ya2  ■  \  -  ^  +  c0  -  \  •  C0Q  ■  c20  +  c30  (37) 


with  the  general  solution 


V„(0)  =  A„  ■  c0  +  • 0 j  •  50  +  y a2  ■  f  c0o  ■  520  +  -  •  c30 

(38) 

Upon  imposing  the  boundary  conditions  stipulated  by  Eqs.  (12), 
the  set  of  constants  in  the  complementary  solution  to 

Eq.  (37)  is  obtained  as 


\B„ 


sin(0/-0o) 


ycr 

v0  -J0/- Y 


j0o  ■  50/  ■  (0/-  0O  —  2c0o  ■  (50/-  50o))  +  — (s0  f  c30o  -  50o  ■  c30/) 


-  V0 '  c9/_ 


y  a 


c0o  •  sOf  (2 sQf  c0o  -  0/)  -  c0/>  50o  ■  (j20o  -  0O)  +  -(c0o  ■  c30/—  c0/  ■  c30o) 


(39) 


Thereafter,  Eq.  (16)  provides  an  analytical  expression  in  closed  form  for  the  total  hydraulic  pressure  field  on  the  upstream  gate  skin,  as 
U— 


,  N  1 dVjfi)  1  s  ya  (  3 

pla,(Q)  = - =  -(-  Aoo  ■  50  +  Ba  ■  c0)  +  —  •  —  0  •  c0  -  s0  +  2c0o  ■  s20  -  —  ■  j30 

a  dQ  a  2  \  4 

Similarly,  when  the  approach  velocity  U  is  very  small  (U  — >0),  Eq.  (33)  for  the  shear  force  function  V(0)  reduces  to 

-y  a2  \  h  1  1  1 

V"  +  V=  —  ■  3-  •  s20  -  -c0  +  -c30 
2  I  a  2  2 


(40) 


(41) 


with  the  general  solution 


,  ,  ya~  I  h  1 

Vo(0)  =  Aa  ■  c0  +  Ba  ■  s0  - - •  0  •  sQ  +  4-  •  j20  +  -  ■  c30 

a  4 


(42) 


Upon  imposing  the  boundary  conditions  stipulated  by  Eqs.  (12),  the  set  of  constants  {A0,i?0}  in  the  complementary  solution  to  Eq.  (41) 
is  obtained  as 


I.  B0  )  sin(0f-0o) 


V0  '  >50/  + 


ya 


ycr 
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h  1  \  hi 
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(43) 


Thereafter,  Eq.  (16)  provides  an  analytical  expression  in  closed  form  for  the  total  hydraulic  pressure  field  on  the  upstream  gate  skin,  as 
U->  0 


pM  =  ~ 


1  dVo(Q)  I  ^  ^  ^  yal  /  h  _  3 


=  -(-  Aa  ■  50  +  Ba  ■  c0)  -  — —  •  0  ■  c0  +  50  +  8~  ■  c20  -  —  •  530 
a  dQ  a  8  \  a  4 


(44) 


Notation 

The  following  symbols  are  used  in  this  paper: 

A;  =  tributary  area  for  the  ith  pressure-cell 
location; 

Aj  =  auxiliary  area  (shaded  in  Fig.  4); 

Ab,/ 1d  =  cross  sectional  areas  of  control  volume  at 
Locations  B  and  D,  respectively; 

A,A',A",A"',AIV 

=  demarcating  points  in  longitudinal  section  of 
control  volume; 

{A,S}  =  set  of  two  integration  constants; 

{A0,B0}  =  set  of  two  integration  constants  when 
approach  velocity  is  very  small; 

{A,,,,^}  =  set  of  two  integration  constants  when 
approach  velocity  is  very  large; 


a  =  gate  radius; 

B ,  B'  =  demarcating  points  in  the  longitudinal  section 
of  control  volume; 

C  =  center  of  gate-skin  circular  arc; 

CD  =  equivalent  viscous  damping  coefficient  for 
drag  force  on  gate; 

Ca  =  equivalent  viscous  damping  coefficient  for 
horizontal  actuator; 

Cc  =  coefficient  of  contraction; 

C*  =  generalized  overall  viscous  damping 

coefficient  for  a  single-degree-of-freedom 
model; 

Cq  =  generalized  viscous  damping  coefficient  for 
internal  energy  dissipation  mechanism  of 
gate  structure  combined  with  effect  of  drag 
force; 
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end 
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M, 

M; 

M. 

M0 


m,  n 


N  or  N(Q) 


N0 


cos(n-Q),  n=  1, ...  ,4; 

demarcating  points  in  longitudinal  section  of 
control  volume; 
differential  operator; 

distance  from  line  of  action  of  resultant  force 
for  pressure  pt  on  gate  skin  to  tube  axis  ( i 

=  1 . 9); 

Young’s  modulus  of  elasticity  for  brass; 
mathematical  expectation  operator; 
base  of  natural  logarithms; 
general  point  on  gate-skin  circular  arc; 
frequency; 

natural  frequency  of  vibration  of  dry-gate 

system  (Hz); 

gate  center  of  gravity; 

half-range  spectral  density  function  for 

pressure  phi=l, ...  ,9\ 

acceleration  of  gravity; 

upstream  head  at  entrance  to  the  control 

volume; 

complex  frequency  response  function  between 

pressure  p,  and  displacement  iq; 

water  depth  of  nappe  at  crest  of  a  sectional 

cantilever  of  gate; 

hydrostatic  critical  depth; 

brink  depth; 

cross-sectional  moment  of  inertia  of  operating 
arms  about  a  horizontal  neutral  axis  through 
their  centroids; 

mass  moment  of  inertia  of  gate  about  tube 
axis; 

imaginary  unit  [\(—  1 )]; 

end  sections  of  operating  arms  adjacent  to 

torque  tube; 

defined  variable  for  3EI/L; 

system  stiffness  matrix; 

axial  stiffness  of  horizontal  actuator; 

generalized  stiffness  coefficient  for  a 

single-degree-of-freedom  model; 

defined  variable  for  1/(1  —  1/tt|2); 

label  for  a  circular  arc  on  gate  skin  parallel  to 

left  edge  of  gate  and  close  to  it,  looking  in 

downstream  direction; 

length  of  gate  operating  arms; 

length  of  a  circular  arc  at  intersection  of  gate 

skin  with  a  vertical  plane,  from  gate  crest  to 

tube  axis,  at  gate  base; 

label  for  a  central  circular  arc  on  gate  skin, 

midway  between  arcs  L  and  R  and  parallel 

to  them,  looking  in  the  downstream  direction; 

resultant  hydraulic  torque  exerted  on  gate 

tube; 

resisting  viscous  moment  at  gate  base; 
resisting  inertial  moment  at  gate  base; 
resisting  elastic  moment  at  gate  base; 
equilibrating  internal  bending  moment  on  gate 
at  elevation  0O; 

subscripts  indicating  that  m  or  n  cycles  of 
vibration  have  elapsed,  respectively; 
internal  axial  force  on  a  gate  general  cross 
section  at  elevation  0; 
equilibrating  internal  axial  force  on  gate,  at 
elevation  0O; 


Pi.P2  =  resultant  forces  of  pressures  on  sides  DD' 
and  AAIV  of  control  volume,  respectively; 

pM  =  atmospheric  pressure; 

pb  =  distributed  vertical  reaction  at  base  of  the 
control  volume; 

Pi  =  pressure  acting  at  ith  pressure-cell  location  on 
gate  skin  (/  =  1 , ...  ,9); 

Pn  =  1'N  approach  hydrostatic  pressure; 
p,(0)  =  total  hydrostatic  pressure  field  acting  on 
upstream  gate  skin  at  general  elevation  0; 
pt0(e),  and  p,„(0) 

=  total  hydrostatic  pressure  field  acting  on  the 
upstream  gate  skin  at  general  elevation  0, 
when  approach  velocity  U  is  very  small,  and 
very  large,  respectively; 

Q  =  point  of  application  of  resultant  force  exerted 
by  gate  on  flow; 

q  =  rate  of  discharge  per  unit  width  of  gate; 

R  =  label  for  a  circular  arc  on  gate  skin  parallel  to 
right  edge  of  gate  and  close  to  it,  looking  in 
downstream  direction; 

R\,R2  —  vertical  and  horizontal  components, 

respectively,  of  force  exerted  by  gate  on  flow; 
rx,r2  =  lever  arms  of  forces  R ,  and  R2,  respectively, 
when  evaluating  their  moments  about  tube 
axis; 

S  =  power  spectral  density  function; 

s  =  circumferential  coordinate  on  a  vertical  plane 
and  along  gate  skin,  with  origin  at  gate  crest, 
and  increasing  downward; 

,sn0  =  sin(n-0),  n=l, ...  ,4; 


U 

u 

V  or  V(0) 

v,m 

V0 

yo(0),yoo(e) 


V 


Vi 

W 

W, 

w 

X,  Y 

*3 


X 

y  and  y2 


a 

7 


time; 

velocity  of  flow  entering  control  volume, 
assumed  horizontal  and  uniform; 
velocity  of  flow  leaving  control  volume, 
assumed  horizontal  and  uniform; 
internal  shear  force  on  a  gate  general  cross 
section  at  elevation  0; 
particular  component  of  general  solution 

y(e); 

equilibrating  internal  shear  force  on  gate,  at 
elevation  0O; 

general  solution  for  Y(0)  when  approach 
velocity  U  is  very  small,  and  very  large, 
respectively; 

abscissa  of  center  of  gravity  of  a  displaced 
gate; 

defined  variable  for  dx-  cp; 
gate  weight; 

weight  of  water  in  the  ;th  component  of 
idealized  control  volume  (1=1, 2, 3); 
ordinate  of  center  of  gravity  of  a  displaced 
gate; 

Cartesian  axes  of  reference; 

lever  arm  of  weight  W3  in  calculating  its 

moment  about  tube  axis; 

centroidal  abscissa  of  auxiliary  area; 

vertical  distance  of  line  of  action  of  force  P2 

from  base  of  control  volume  and  from  tube 

axis,  respectively; 

angle  that  gate  operating  arms  make  with 
vertical; 

unit  weight  of  water; 
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7;  = 
A  = 

8  = 
K  = 

8(f)  = 
e  = 

T1  = 
0  = 

0^  and  0O  = 

0[  and  02  = 

M-  = 

vn  and  vn+m 


p  = 
2  = 
(T2  = 
to  = 

tp  and  ip L  = 
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